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Abstract 

We discuss how to extract the spectroscopy of quantum chromodynamics 
(QCD) in the pure gauge sector from the Hamiltonian lattice field theory 
approach. The recently developed truncated eigenvalue equation method is 
applied to the estimation of the scalar glueball ++ and masses in the 
(2+l)-dimensional case. These masses reach the constant values in a scaling 
region as required by the renormalizability. 
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QCD predicts the existence of glueballs, the gluonic bound states formed 
through strong interactions of gluons. Precise calculations of the glueball 
spectrum will also test decisively the validity of QCD and give guide to expe- 
rimentalists in future particle searches. One of the most popular and practical 
techniques for the nonperturbtive determination of the glueball spectroscopy 
is to do numerical simulations on the Euclidean lattice. A lot of progress in 
this direction has been made, giving increasingly accurate estimation of the 
glueball masses. 

Here we would like to discuss an alternative approach: the Hamiltonian 
formulation, which by solving the lattice Schrodinger equation can directly 
provide not only the glueball masses from the eigenvalues, but also the pro- 
files, i.e., the wave functions for the ground state and excited states. 

QCD in 3 dimensions (QCD 3 ) [?, ?] is a nontrivial SU(3) gauge theory, 
which shares a lot of important features of QCD 4 , but is much more simplified 
due to lower dimensionality and superrenormalizability. Furthermore, it is 
hopeful that some techniques and findings in QCD 3 may be generalized to 
its 4 dimensional physical relative. 

In a previous paper [?], we have investigated the long wavelength struc- 
ture of the ground state of QCD 3 , and evaluated its vacuum wave function. 
It has been demonstrated that our lattice vacuum wave function has the 
correct continuum limit and nice scaling behavior for the coefficients in the 
continuum vacuum wave function. This is the first step towards the under- 
standing of the structure of the glueball and hadrons. It is worth mentioning 
that the ground state properties of the lattice SU(2) pure gauge models, 
being closely related to the confinement phenomenon, have been extensively 
studied both numerically [?, ?, ?] and analytically [?, ?, ?, ?, ?, ?]. What 
is new in [?] is the exploratory study of the realistic gauge group SU(3) 
using the recently developed truncated eigenvalue equation method [?, ?] 
and some new prescription scheme for the classification of the graphs in the 
wave function, which lead to rapid convergence to the continuum limit. 

The purpose of this paper is to give a first estimation of the scalar glueball 
masses in QCD 3 as well as their wave functions. 

Concerning the notations, truncated eigenvalue equation method and 
classification scheme for the graphs, we refer to [?] for further details and 
here just give the most relevant ones. On a discrete space and continuous 
time lattice with temporal gauge, the dynamics of the pure gauge system is 
described by 
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H = fj2 EfEf ~— 2 Yl Tr(U p + W p - 2), (1) 
2a J ag z ^ 

which becomes the Yang-Mills Hamiltonian with a renormalized charge e 
in the continuum limit (a — > 0, or equivalently f3 = 6/g 2 — > oo). The 
superrenormalzabilty of the theory in 2+1 dimensions implies that the bare 
coupling g and the lattice spacing a have a simple relation g 2 = e 2 a. 

The ground state can be written in the form of |fi) = exp[R(U)] |0), with 
|0) being the fluxless bare vacuum and R(U) being a linear combination of 
gauge invariant gluonic operators such as the Wilson loops. In [?], we have 
illustrated how to obtain R(U) and the vacuum energy by a truncation 
method developed in [?], where R(U) is expanded in order of graphs R(U) = 
J2iRi(U) } and Ri(U) is classified according to the scheme described in [?], 
with i being the order of the graphs. Up to the third order, they are 

R^U) = C X G X + h.c, 

8 = 1 
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R3(U) = J2C3,G 3 , t + h.c. } 



1=1 



which graphs G\ } G^,i and Gz,i are given in Fig. 1 (see also [?]), and the co- 
efficients Ci, C2,i and Cz,i are determined by solving the eigenvalue equation 
at the corresponding order n: 

n c\ 

E [E h R t (U)][E h R 3 (U)]}--J2Tr(U p + Ui) 

l i i+j<n ^ V 



2a , , 

= -en- 3 

9 

To calculate the glueball masses, an important procedure is to construct 
appropriate states for the particles. The wave functions of these states are 
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created by gauge invariant operators with the same quantum numbers as 
these glueballs acting on the vacuum state. Here for simplicity we concentrate 
on two scalar glueballs with quantum numbers J pc = ++ and J pc = . 
The wave function of a glueball is of the form 

which is so chosen that it is orthogonal to the vacuum state |fi). The creation 
operator of the state Fill) is a linear combination of a complete set of graphs 
with the given quantum number of the glueball. The gap or the mass of the 
glueball Mjpc = Ae = ep — €q can be obtained by solving the Schrodinger 
equation 

H\F) = e F \F), (5) 
which is reduced to the eigenvalue equation 

^{[Et, [E h F(U)\] + 2[E h F(U)][E h R(U)]} = ^F(U). (6) 

Again, Fill) is expanded in order of graphs as 

F(U) = J2m), (7) 

i 

and (??) is truncated to the nth order as 

n cy a n 

EP'.[^E f .( f/ )]] + 2 E [E h F i {U)][E h R j {U)]} = —Y t F{U). (8) 

/ i i-\-j<~n $ i 

We choose the operators for the glueball ++ as 
Ff\U)=f^G 1 + h.c, 

*f + (U)=j2f£ + G2,i + h.c, 

2 = 1 

29 

Ft +(U) = J2f£ + G^ + h.c, (9) 

8 = 1 
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which are the same set of graphs G\ } G^,i and Gz,i as those in (??) but with 
different coefficients. Substituting (??) into (??), we obtain 36 nonlinear 
equations for /° ++ , f° \ + ■> an d f% ; + ? which can be solved numerically. 
The operators for the glueball are 

Fr(u) = frGi-h.c. } 

^2 (U) = G 2 ,l + /°,3 ^2,3 + /°,4 ^2,4 ~ h.C, 

^3 "(^) = E/3 rG'°7-/ i .C., (fO) 

where the coefficients are similarly determined and G® i stand for the ope- 
rators with J pc = chosen from G^^. 

For large f3 = 6/(/ 2 , there should be a so called scaling region where the 
physical quantities become constants. Equivalently, by dimensional analysis, 
in the lattice unit the inverse correlation length or the dimensionless masses 
aMjpc should be proportional to g 2 } i.e., 

aMjpc Mjpc 

5 > 5— = const., (11) 

as required by the renormalizability of the theory. 

The results for aM ++ / g 2 and aM — / g 2 are shown in Fig. 2. As one 
sees, there exists indeed a scaling window f3 £ [5, 8) for ++ and f3 £ (6, f 2] 
for where we can extract their masses (indicated by the dot lines in Fig. 
2). From the data on the dot lines, we get 

it: K h6m (12) 

To summarize, using the techniques and information on the vacuum wave 
function in [?], and solving the lattice Schrodinger equation truncated to the 
third order, we have obtained the glueball wave functions and masses in 
QCD 3 . The results for the mass spectrum show sensible scaling behavior 
in some intermediate or weak coupling regions. (In comparison, the QCD 4 
glueball masses from the Monte Carlo simulations in the literature are ex- 
tracted within f3 £ [5.9,6.2].) Of course, when approaching the continuum 
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limit a — > or equivalently f3 — > oo, as the correlation lengths in the lattice 
units will be divergent, a conventional wisdom is to include higher and higher 
orders of graphs to better represent the glueballs. This is quite similar to the 
situation in the Monte Carlo simulation on the Euclidean lattice where larger 
and larger lattices are required in the weak coupling regime. It may also be 
possible to improve the scaling using a better graphic classification (mixing) 
scheme. This issue and the higher order effects on the physical observables 
will be further investigated. 
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Figure Captions 
Fig. 1.1. First order graph in RiU). 
Fig. 1.2. Second order graphs in RiU). 
Fig. 1.3(a). Third order graphs in RiU). 

Fig. 1.3(b). Third order graphs (continuation of Fig. 1.3(a)) in RiU). 

Fig. 2. aM ++ 1 g 2 and aM — / g 2 as a function of f3. The dot lines indicate 
the corresponding mean values in the scaling region. 
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Fig. 2 
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Fig. 1.3(a) 
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Fig. 1.3(b) 



